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1. Introduction 

The pivotal role played by gauge field theories in a unified description of fundamental 
interactions proposed one of the most challenging questions of modern high energy 
theoretical physics: if Nature likes so much gauge symmetry why gravity cannot fit into 
such an elegant and, would be, universal blueprint? 

Before the advent of string theory this was a question without an answer. After that, 
it became clear that all field theories, including Yang-Mills type models, must be seen 
as low energy, effective approximations of some more fundamental theory where the 
dynamical degrees of freedom are carried by relativistic extended objects. Furthermore, 
even the low-energy effective gauge theories require a "stringy" approach in the strong 
coupling regime, where standard perturbation series breaks down. Color confinement 
in QCD is a remarkable example of a phenomenon where the string tenet meets gauge 
symmetry. The stringy aspects of confinement have long been studied, but are not yet 
fully understood^. Several different models have been proposed as phenomenological 
descriptions of the quark-gluon bound states, including color flux tubes, three-string of 
various shapes, bag-models [|J . To promote some of them to a deeper status one would 
like to derive extended objects as non-perturbative excitations of an underlying gauge 
theory [[J. The first remarkable achievement of this program was to obtain stringy 
objects from SU(N) Yang-Mills theory in the large-iV limit ||. 

These results have been extended to the case of a self-dual membrane in the SU (oo) 
Toda model 0. In a nutshell, the problem is to establish a formal correspondence be- 
tween a Yang-Mills connection, A l j^(x), and the string coordinates, X M (cr° , a 1 ), i.e. 
one has to "get rid of" the internal, non-Abelian, indices i, j and replace the spacetime 
coordinates x M with two continuous coordinates (a , cr 1 ). With hindsight, the "recipe" 
to turn a non-Abelian gauge field into a set of Abelian functions describing the embed- 
ding of the string world-sheet into target spacetime can be summarized as follows: 

a) transform the original field theory into a sort of "matrix quantum mechanics" ; 

b) use the Wigner-Weyl-Moyal map to build up the symbol associated to the above 
matrix model; 

c) take the " classical limif of the theory obtained in b). 

Stage a) requires two sub-steps: 
al) take the large- N limit, i.e. let the row and column labels i, j to range over arbitrarily 
large values; a2) dispose of the spacetime coordinate dependence through the so called 
"quenching approximation", i.e. a technical manipulation which is formally equivalent 
to collapse the whole spacetime to a single point. 

After stage a) the original gauge theory is transformed into a quantum mechanical 
model where the physical degrees of freedom are carried by large coordinate independent 
matrices. Then, stage b) associates to each of such big matrices its corresponding 

+ There are many interesting reviews on this problem, as (lj, S, (3) : ... and we apologize for omitting 
many other good ones. 
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symbol, i.e. a function defined over an appropriate non-commutative phase space. The 
resulting theory is a deformation of an ordinary field theory, where the ordinary product 
between functions is replaced by a non- commutative * -product. The deformation 
parameter, measuring the amount of non-commutativity, results to be 1/N, and the 
classical limit corresponds exactly to the large- N limit. The final result, obtained at 
the stage c), is a string action of the Schild type, which is invariant under area-preserving 
reparametrization of the world-sheet. 

More recently, we have also shown that bag-like objects fit the large-iV spectrum of 
Yang-Mills type theories as well, both in four || and higher dimensions. We started 
from the Yang-Mills action for an SU( N ) gauge theory supplemented by a topological 
term 

S = -Jd*x( J*- TrF^ + 6 -^r L Tri ^ ) C 1 ) 

and went trough the steps from a) to c). As a final result we obtained the following 
action 

w = ~Yq J*<t{X",X v } pb {X ii ,X v } pb + 

-«W, / £sX» {X v ,X" ,X a } NPB , (2) 
Jan 

where 

{X» ,X"} PB = e mn d m X»d n X» 

is the Poisson Bracket and 

{X v ,X" ,X°} NPB = ^ k d i X v d j X"d k X <r 

is the Nambu-Poisson Bracket. The first term in (0) describes a bulk three-brane, or 
bag, which in four dimensions is a pure volume term characterized by a pressure /j,q. 
All the dynamical degrees of freedom are carried by the second term in (0), where k is 
the membrane tension; this term encodes the dynamics of the boundary, Chern-Simons 
membrane, enclosing the bag. Tracing back the bulk and boundary terms in the original 
action ([[]) it is possible to establish the following formal correspondence 

"glue" « — > bulk 3-brane 

"instantons" « — > Chern-Simons boundary 2-brane. 

This scheme, which has been generalized to Yang-Mills theories in higher dimensional 
spacetime 0, points out that not only strings but bag-like objects fit the large- iV 
spectrum of SU(N) gauge theories. However, a non-trivial dynamics for these spacetime 
filling objects comes only from boundary effects, described by Chern-Simons terms in 
the original gauge action. 

Against this background, we would like to investigate the existence of non-topological 
membrane-like excitations in the large large- iV spectrum of a SU(N) Yang-Mills theory 
in four dimensions. Clues suggesting the existence of these objects come from earlier 
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Abelian models |TD) and the recent conjectures about M-Theory. M(atrix) theory is 
the, alleged, ultimate non-perturbative formulation of string theory. More in detail, two 
models have been constructed as possible non-perturbative realization of Type IIA [] I 



and Type IIB |T2[ string theory. The matrix formulation of Type IIB strings is provided 



by a large- N, 10-dimensional super Yang-Mills theory reduced to a single point 

S IKKT = -jTr[A^A u ] 2 + .... (3) 

The dots refers to the fermionic part of the action which is not relevant to our discussion. 
The model (|3|) has a rich spectrum of extended objects. Our investigations in || and 
H has been initially triggered by the formal analogy between in 10D and quenched 
Yang-Mills theory in AD. 

Matrix description of Type IIA strings is given in terms of 0-branes quantum mechanics 

-\[*. (4) 

where i = 1, . . . , 9. Again, the 0-branes matrix coordinates can be seen as Yang-Mills 
fields reduced to a line. 

In this paper we would like to "reverse the path" leading Type IIA strings to the matrix 
model (U) and show how a 3-D version of (§) can be obtained from the canonical for- 
mulation of an SU( N ) Yang-Mills theory through a modified quenching prescription. 
Then, we are going to extract a non-relativistic, dynamical 2-brane from the large- N 
spectrum of the model by following the procedure introduced in and |§. 

The paper is organized as follows: in section |2] we start from SU(N) Yang- 
Mills theory in four dimensions and obtain a corresponding matrix theory through the 



quenching approximation; two different type of quenching are discussed in section 2.1 



and section |2.2| ; in section g we study the large- N limit of the matrix model introduced 
in section ^]2| and use the Weyl-Wigner-Moyal map to get the action for a membrane; 
we conclude the paper by computing the mass spectrum of this membrane in the WKB 
approximation. 

2. Yang-Mills Theory as a Matrix model 

In the introduction we referred to the supposed relation between confinement and 
extended excitations of the Yang-Mills field. Recently, an even deeper and more 
fundamental relation between branes and Yang-Mills fields has been conjectured in 
the framework of M-theory [13| . Non-perturbative formulations of string theory require 



the introduction of higher dimensional solitonic objects satisfying Dirichlet boundary 



conditions |14]]. Dirichlet-branes are formally described by non-commuting matrix 
coordinates. Thus, non-perturbative string theory, or M-theory, is conveniently written 
in terms of matrix dynamical variables. The corresponding low energy, effective, 
supersymmetric Yang-Mills theory is derived through an appropriate compactification 
procedure of the original matrix model [11], Jl5|| . 
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In this section we are going to follow a similar path, but in the opposite direction: we 
shall start from a SU(N) gauge theory in four dimensional spacetime and build a matrix 
model. Our final purpose is to show that the spectrum of such a QCD matrix model 
contains dynamical membrane type objects. We thus start from the Yang-Mills action 

S YM = [dt [ d 3 xC YM (F ,A) (5) 
J Jv H 

defined in terms of the Lagrangian density 

,-^Tr {F, u f - -1 
The covariant derivative is defined as usual 



C YM (F,A) = —Tr (F^f - —TrF^D^ . (6) 



DfrAv] = dfrAv] + % [An, A v ] , (7) 
but the SU (N) Yang-Mills Lagrangian £ym {E, A) is written in the first order 



formulation: thus 



A — Aarpa jp - no rpa 



are independent vector and tensor fields respectively^ valued in the Lie algebra defined 
by the commutation relations 

The integration volume Vjj will be specified later on. The form (0) is appropriate for an 
Hamiltonian formulation of the action, as it is required by the quenching approximation 
that we shall apply in the next section. 

Starting from the Lagrangian density we can, as usual, introduce the canonical 
momentum and Hamiltonian 

, _ <9£ Y M 1 jptm (q\ 

E = nn^T - — oT~2 W 



dd t Ai 2g t 



o 



H = Tr (E'dtAi) - C (9) 
and rewrite (||), in terms of phase space variables, as: 

Cym= ~ TT^Tr (F tm ) 2 - -^TrF' m («9 t A, + d t A t - i [A t: A t ]) 



J^_Tr (F tm ) 2 

+ -^Tr (F mn ) 2 - ^TrF mn D [m A nu] 
4g 2 2g 2 



~ f^Tr (E'f + i T r (E^A,) - ^Tr (A tDl E) - ^Tr (D [m A n] ) 2 . (10) 
Accordingly, the phase space action reads 



S= [dt [d 3 x [E%Ai - Ho] (11) 
J Jv H 

= Jdt Jd 3 x ^Tr (Ef- ^Tr (E%A t ) + ^Tr (A t A#) + ^Tr (D [m A n] ) 



The metric signature is ( h + + ) and in our notation matrices are denoted by boldface letters. 
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Let us remark that A t enters linearly in the canonical form of the action (|T2| ) and plays 
the role of Lagrange multiplier enforcing the (non-Abelian) Gauss Law: 

5S 



5 A t 







DiE 1 = 0. 



(13) 



Thus, solving the classical field equation for A t is equivalent as requiring E l to be 
covariantly divergence free in vacuum. Thus, inserting the solution of the Gauss 
constraint (|IB|), the action for E l becomes 

2 - 1„ . s N 



S- 



YM 



9o 
2N 



dt / d 3 x 
Jv H 

We will now go on with the quenching procedure. 



Tr(^) --Tr(Ed t A t ) + —Tr(D [m A n] y 



(14) 



2.1. Spacetime Quenching — > IKKT-type model 

In order to provide the reader a self-contained derivation of our model, let us briefly 
review how the quenching approximation works in a simple toy model pi. Consider 



the following two-dimensional model of matrix non-relativistic quantum mechanics 



S 



d 2 x Tr 



1 9 1 

- (d M) 2 -- 
2 y 1 2 



d\M ) 2 -V(M) 



(15) 



where M( x° , x 1 ) is an Hermitian, 2x2 matrix and V ( M ) is an appropriate potential 
term and we suppose that the system described by M is enclosed in a (one-dimensional) 
"spatial box" of size a. Quenching is an approximation borrowed from the theory of 
spin glasses where only "slow" momentum modes are kept to compute the spectrum of 
the model. Slow modes are described by the eigenvalues of the linear momentum matrix 
P, while the off-diagonal, "fast" modes can be thought as being integrated out. Thus, 
M( x ) can be written as 

M(x° ,x x ) =exp (iPx 1 ) M(x°)exp(-iPx 1 ) (16) 

and the spatial derivative d\M{x° ,x x ) becomes the the commutator of P and M 



diM = i [P ,M 
so that the action becomes 

S = a I dx° Tr 



- («9 M) 2 + - [P ,M} 2 -V(M 



(17) 



(18) 



in which the the x 1 dependence of the original matrix field has been removed. 
Quenching can be applied to a Yang-Mills gauge theory by taking into account that in 
the large- iV limit SU(N) — > U (N) and the group of spacetime translations fits into the 
diagonal part of U(oo). By neglecting off-diagonal components, spacetime dependent 
dynamical variables can be shifted to the origin by means of a translation operator 
U(x): since the translation group is Abelian one can choose the matrix U(x) to be a 
plane wave diagonal matrix |T 



U ab (x) = 5 ab exp (iq a 



I. > L 



(19) 
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where q a M are the eigenvalues of the four-momentum q^. Then 

A„(x) = exp (-iqr^) A„(0) exp (iq^) = U\x)AfU{x) 

and in view of the equality 

D„A V = iU\x) [q„ + A®, A v ] U(x), 

which when antisymmetrized yields 

D { ,A V] = iU\x) [q, + Af,q v + Af] U(x) = iU\x) [A®, A$>] U(x), 

we can see that the translation is compatible with the covariant differentiation, so that 

F^x) = exp {-iq^) F^(0) exp (iq^) = U\x)F$U{x). 

Once the original gauge field theory is turned into a constant matrix model, we still need 
to dispose of the spacetime volume integration. The gluon field is spatially confined 
inside a volume Vh comparable with the typical size of an hadron. Thus, for any finite 
time interval T we can replace the four-volume integral by 



T 



dt I d 3 x — ► TV„ 



H 

JV H 



and the quenched action becomes 

4^.=^^^^ , (20) 

which is the first order formulation of the IKKT-type action in four spacetime 
dimensions [12]]. The usual second order formulation is readily obtained by solving 
for Fj^) in terms of AjL 

FV>=i[AU>,AF] 
and substituting back this result into 



4 q i- red . ^ S%* T = /3V h JLtt ([A^,A^]) 2 . (21) 

The string-like excitations of this model and the relation between large-iV gauge 
symmetry and area-preserving diffeomorphism have been investigated in several papers 
||. More recently, we found that not only strings are present in the large- N spectrum 
of (|21~D but also spacetime filling, bag-like objects [§], for which a non-trivial boundary 
dynamics was found through the addition of topological terms to the original Yang-Mills 
action. Here, we would like to explore a different route leading in a more straightforward 
way to a dynamical brane action. From this purpose we need to introduce a different 
quenching approximation, which we discuss in the next section. 
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2.2. Spatial Quenching — > BFSS-type model 

Instead of shifting A M (i, x) to a single point, as we did in the final part of the previous 
section, we translate the matrix gauge field to a fixed time slice by means of the conserved 
three-momentum q. As in the previously discussed case 

Ai(t,x) = exp (-iqiX 1 ) Ai(t)exp (iqix 1 ) = U\x) Ai(t)V '(x) 
A t (t,x) = exp (-iqiX*) A t (t) exp fax') = U\x)A t {t)U{x), 

the translation operation commutes with the covariant differentiation since 



DiA r 



iU\x)[ qi + Mt),A n ]U(x) 



implies 



iU\x) [q, + Ai(t), q n + A v {t)} U(x) 
and for the conjugate momentum we also get 



U\x)E,{t)U{x) = iU\x) lA^it) ,A$>( t ) 



U(x). 



Enclosing the system in a proper quantization volume Vh while keeping the time 
integration free 

I d 3 x > Vh 

Jv h 



H 



we get 



S = V, 



H 



dt 



2N 



Tr {E\t))' 



1 



Tr E 



idAj 
dt 



N 



rTr 



A^\A^ 



1 \ 2 



which is the action in the first order formulation; by substituting for E l (t) its expression 
in terms of the vector potential we obtain, in the second order formulation, an action 
quite similar to the one for the bosonic sector of the BFSS model describing a system 
of iV .DO-branes in the gauge A = 0: 

2 



s 



BFSS 



N 

V H — I dt 

9o 2 




(22) 



the only difference is the range of the spatial indices: we are working in three rather 
than nine spatial dimensions. An action of the form ( f22] ) can also be obtained from 
monopole condensation and toroidal compactification [p~8|j . 



3. Non-commutative Phase Space 

To match the large-iV SU(N) gauge theory with some appropriate brane model we 
have to bridge the gap between non-commuting Yang-Mills matrices and commuting 
brane coordinates. Since the world-trajectory of ap-brane is the target spacetime image 
= A^(a°, a 1 ,..., a p ) of the world manifold S : a m = (a , a 1 ,... cr p ), A M belonging 
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to the algebra A of C°° functions over E, to realize our program we must deform A to 
a non- commutative "starred" algebra by introducing a ^-product. The general rule is to 
define the new product between two functions as (for a recent review see [fl9| ): 

f*g = fg + hP h (f ,g) , (23) 

where Pn( f ,g) is a bilinear map P% '■ A x A — > A. h is the deformation parameter 
which is often denoted by the same symbol as the Planck constant to stress the analogy 
with quantum mechanics, where classically commuting dynamical variables are replaced 
by non-commuting operators. In our case the role of deformation parameter is played 
by 

For our purposes, we select £ = M 2n and choose the Moyal product as the deformed 
*-product 

d 2 



f(a)*g(a) 



cxp 



.h , 

l— OJ 

2 



f(cr)g(0 



(25) 



£=<7 



where u mn is a non-degenerate, antisymmetric matrix, which can be locally written as 

In. X n 



-I 



(26) 



The Moyal product ( p5|) takes a simple looking form in Fourier space 

F{°) * = | ^ exp (^ mn a m r) F (| + C) G (| - C) , (27) 

where F and G are the Fourier transform of / and g. Let us consider the Heisenberg 
algebra 

[K rn , P 1 } = iU ml ; (28) 

Weyl suggested, many years ago, how an operator Op(K, P) can be written as a sum 
of algebra elements as 

1 



Or 



(27T) 



d n pd n kF(p, k) exp (ip m K m + i^P 1 ) 



(29) 



The Weyl map (|29|) can be inverted to associate functions, or more exactly symbols, to 
operators 



F(q ,k 



- — — exp ( — ikt ) ( q + n- 



O f (K ,P 



q — h- 



(30) 



moreover it translates the commutator between two operators t/, V into the Moyal 
Bracket between their corresponding symbols U(o~), V(a) 



ill 



[U, V] <— > {W, V} MB = - (U * V - V * W) 



and the quantum mechanical trace into an integral over Fourier space 

(2vr) n Tr w O F (K,P) i — > J d n pd n k F(q ,k) = J d 2n aF( 



a 



(31) 
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A concise but pedagogical introduction to the deformed differential calculus and its 



application to the theory of integrable system can be found in [20 



We are now ready to formulate the alleged relationship between the quenched model 
(f22|) and membrane model: the symbol of the matrix Aj is proportional to the 2n- 
brane coordinate X^a 1 , . . . , cr 2n ). Going through the steps discussed above the action 
gBFSS transforms into its symbol jy BFSS : 



1 dAi dA % {Ai,A n } MB *{A\A n } MB 
* —rr + 



(32) 



2 dt dt 4 

The action ([32]) is manifestly Lorentz non-covariant, as it is expected (the covariant, 
super symmetric, higher dimensional version of the action (^) is discussed in ||21|| ). 
The adopted quenching scheme explicitly breaks the equivalence between spacelike and 
timelike coordinates. Accordingly, our final result takes a typical "non-relativistic" look. 
Up to now we have not fixed the Fourier space dimension n. To give A 1 the meaning of 
embedding function, we have to choose 2n < D — 1, where D is the target spacetime 
dimension. To match QCD in four spacetime dimensions we set n = 1. In this case 

u mn = e mn =(^! I J ( 33 ) 
and we rescale the Yang-Mills charge and field| as 

4 ~ 4- (34) 
9o 2 9ym 

A' i — ► V~ 2IZ X { . (35) 

Since the "glue" is supposed to be confined inside an hadronic size volume Vh, we can 
assign to gyu the standard value at the confinement scale 

^XM~0.18. (36) 
An 

Finally, if A^ 3> 1 the Moyal bracket can be approximated by the Poisson bracket 

W^U— ^ {**>*'}PB 
jj For the sake of clarity, let us summarize the canonical dimensions in natural units of various 
quantities: 

[A^(x)] = [AW] = (length)" 1 
[F^ a (x)} = [F^W] = (length)- 2 

[a m ] = (length) = 1 , [t] = length , [/3] = length 
[go] = [qym] = (length) = 1 
[V H ] = (length) 3 
[X 1 ] = length 

[ Mo ] = (length)- 1 , [a] = (length)" 5 . 



The QCD Membrane 



11 



and (1321) takes the form Wl 



^nr 2 = J dtd 2 a 
where /xo and a are defined by 



1 dX fc dX fc a 
2^° (it dt + 4 



+ T {^,X J } PB {X\^} 



PB 



(37) 



-1/3 



i3 



-5/3 



a: 



H 



2vr^ M " 2tt^ m 



Moreover from the definition of the Poisson bracket 



{X\X j } 



PB 



we can compute 



{X 1 ,Xi} pB {X { ,X 3 } pB = 2-det(d m X k d n X k )=2 1 



so that the action 



q p=2 



can be rewritten as 



dtd a 



1 dX k dX k a , ^ , _ ^ fe „ _ _ \ 

■A^o-jr -JT- + 2 det ( dm x dn Xk ) 



(3f 



_2 r " u dt dt 

The first term in (|38D is a straightforward generalization of the kinetic energy of a 
non-relativistic particle; the second term represents the "potential energy" associated 
to the elastic deformations of the membrane. The action ( [38|) still displays a residual 
symmetry under area preserving diffeomorphisms, leaving only one dynamical degree 
of freedom describing transverse oscillations of the membrane surface|tt|. If the action 
(|38|) has any chance to provide a membrane model of hadronic objects, then it must be 
able to provide at least the correct order of magnitude of hadronic masses. Our model 
does not take into account spin effects, therefore it is consistent to look for spherically 
symmetric configurations. Again this is a sort of quenching even if of a more geometric 
type. Infinite vibration modes of the brane, corresponding to local shape deformations, 
are frozen and the dynamics is reduced to the "radial" breathing mode alone. This 
kind of approximation, commonly called "minisuperspace" approximation, is currently 
adopted in Quantum Cosmology, where it amounts, in practice, to quantize a single scale 
factor (thereby selecting a class of cosmological models, for instance, the Friedman- 
Robert son- Walker spacetimes) while neglecting the quantum fluctuations of the full 
metric. The effect is to turn the exact, but intractable, Wheeler-DeWitt functional 
equation [23| into an ordinary quantum mechanical wave equation |24]. As a matter of 
fact, the various forms of the "wave function of the universe" that attempt to describe 
the quantum birth of the cosmos are obtained through this kind of approximation 



or modern refinements of it |26] . This non-standard approximation scheme was applied 



ff We have assumed that X 1 are three spacelike coordinates. By relaxing this assumptions we can give 
( |38| ) a slightly different physical interpretation. If all three X 1 are considered as transverse directions, 
then ([38]) can be seen as the light-cone gauge action for a bosonic brane in a 5-dimensional target 
spacetime. 
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to a relativistic membrane in the seminal paper by Collins and Tucker p7|, and since 



then it has been used several times [28], including the case of self-gravitating objects 



Following p?fl , we parametrize the membrane coordinates as follows 

X 1 = R(t) sin 9 cos 
X 2 = R(t) sin 6 sin 
X 3 = R(t) cos 9 



(39) 



and the transverse, dynamical degree of freedom corresponds to R. The metric ^ a b 
induced on the membrane by the embedding (|3~9"1) is: 



lab = diag (R 2 (t), R 2 (t) sin 2 9) , det ( 7mn ) = R\t) sin 2 9 
The corresponding action turns out to be 



S 



Ldt = vr 2 dt 



1 / dR\ a 4 
2^Vdi) + ~A R 



accordingly the momentum conjugated to the only dynamical degree of freedom is 

Pr = 



dL 2 dR 

d(dR/dl) ~ 71 



from which the Hamiltonian can be calculated as 



K dt 



Then the action in Hamiltonian form is 

dR 



1 P 2 + ^R 4 



2vr 2 /i 



S 



dt 



a 



27T 2 /i H 4 



The above results allow one to compute the hadronic mass spectrum from the spherical 
membrane Schrodinger equation 



d 2 



+ 



an 



R 4 



2n 2 [i dR 2 ' 4 
with the following boundary conditions 

#(0) = 
lim V(R) = 0. 



q(R) = M n ^(R) 



(40) 

(41) 
(42) 



The lowest mass eigenvalues can be evaluated numerically [|30| or through WKB 
approximation 



(27rV M n ) 



1/2 



47T 



/?(l/4,3/2) \8V2V H g 2 



'YM 



1/3 



n 



2/3 



(43) 
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where, (3 is the Euler /^-function: £(1/4,3/2) = r(1 ffJg /2) . The WKB formula 

2/q 1/3 

gives a mass scale of the correct order of mag nitude, M n oc ^ g^ V H - 1GeV - More 
sophisticated estimates of the glueball mass spectrum, including topological corrections 
are not very different from the values given by ([Ef). Thus, we conclude that the 
QCD membrane action fl38|) encodes, at least the dominant contribution, to the gluon 
bound states spectrum. Hopefully, an improvement of this result will come from an 



extension of the minisuperspace approximation along the line discussed in |]3"2" |, where a 
new form of the p-brane propagator has been obtained. 
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